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- - , Abstract 

m . 

■ We generalize the Harish-Chandra-Itzykson-Zuber and certain other integrals (Gross- 

i Witten integral and integrals over complex matrices) using the notion of tau function of 

^ ' matrix argument. In this case one can reduce multi-matrix integrals to the integral over 

Q^' eigenvalues, which in turn is a certain tau function. We also consider a generalization of 

, the Kontsevich integral. 

m ■ 

: 1 Introduction 

^ ■ The study of different models of random matrices is one of the most interesting topic in math- 
! ematical physics (see [1],[2],[3],[4] for a review). Models of random matrices are applied to 
, various topics starting with quantum field theory and condensed matter physics ending with 
combinatorics, algebraic topology and number theory. The purpose of the paper is to generalize 
' some known solvable matrix integrals. By the solvable matrix integrals we mean those which 
^ ■ can be presented as tau functions and which can be evaluated by the method of orthogonal 
\Q , polynomials. Let us mark that in papers [5], [6] by solvability of matrix model the author means 
O . a weaker condition - a model considered to be solvable if it admits a reduction of number of 
^ I integrations (from the order ~A^^ to the order A^) via a character expansion method without 



O 



X 



CN any relation to the tau function topic. 

Today it is not surprising that the notion of tau function, introduced by Sato school (see a list 
of references in [7]), is widely used in statistical physics and quantum field theory as a partition 
^ ! function of different models (as an example see the review [8]). The special class of tau functions, 
■ which we shall exploit in the present paper and which we call tau functions of hypergeometric 
type, is not an exclusion. Tau functions of hypergeometric type are parameterized by a function 
r on the lattice and are denoted by r^. (The name hypergeometric tau function [9] is due to the 
fact that the simplest specialization, r is rational, yields the hypergeometric function of matrix 
argument [10], [11].) This is a type of tau functions one may meet in many different problems, 
and special examples of these tau functions one can find in [12], [13], [14], [15], all these examples 
one may obtain by a specialization of r. The most of well known matrix integrals actually belong 
to this class of tau-functions [16], [17], [18]. (To be precise we deal with asymptotic expansion 
of matrix integrals and tau functions). 

In the present paper we show that the integrals of tau functions of hypergeometric type of 
matrix argument are again tau functions of hypergeometric type. We notice that all known 
matrix integrals such as one-matrix model [19], Itzykson-Zuber integral [20], [21], two-matrix 
matrix model [21], Gross- Witten one plaquette model [22] may be considered as integrals of the 
simplest tau function (so-called vacuum tau function). The determinant representation of the 
tau functions of hypergeometric type allows to obtain generalizations of solvable multi-matrix 
models and a certain generalization of Kontsevich integral. 
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In the last part of this introduction we review some things from sohton theory. 
Soliton theory. KP hierarchy of integrable equations [23], [7], which is the most popular 
example, consists of semi-infinite set of nonlinear partial differential equations 

9tm'^ = K^M , m = 1, 2, . . . , (1) 

which are commuting fiows: [9^^, d^^] u — 0. The first nontrivial one is Kadomtsev-Petviashvili 
equation 

d,,u = ^dlu + Idi^'dlu + ^dt^ , (2) 

which originally served in plasma physics [23] , now plays a very important role both, in physics 
and in mathematics. Another very important equation is the equation of two-dimensional Toda 
lattice (TL) carefully studied in [24] 

ai^,9q0n = e'^"-i-'^"-e'^"~'^"+i (3) 

This equation gives rise to TL hierarchy which contains derivatives with respect to the higher 
times ti,t2, ■ ■ ■ and tl, t^, . . . . 

The key point of the soliton theory is the notion of tau function, introduced by Sato (for 
KP tau-function see [7]). The tau function is a sort of a potential which gives rise both to TL 
hierarchy and KP hierarchy. It depends on two semi-infinite sets of higher times ti,t2, ■ ■ ■ and 
tl, . . ., and on discrete variable n: r — T{n, t, t*). More explicitly we have [7], [24]: 

u = 2dl log r(n, t, t*) , 0„(t, t*) = - log ^^""^^Y'^P (4) 

In the soliton theory so-called Hirota-Miwa variables x, y, which are related to the higher 
times as 

mtm^Y.< ^ mC^J^vT, (5) 

i i 

are well-known. Any tau function is a symmetric function in Hirota-Miwa variables, the higher 
times mtm, are so-called [25] power sums. 

It is known fact that a typical tau function may be presented in the form of double series 
over partitions in Schur functions [26]: 

r(n,t,r) = ^i^A^SA(t)s^(t*) , (6) 

where the coefficients Kxfj. solve special bilinear equations [7]. 

T functions of hypergeometric type. Let us consider a function r which depends on a 
single variable n, the n is an integer. We suppose r{n) to be finite for integer n. Given partition 
A, we define 

'^>^{^) = n r{x + j-i) (7) 

Namely, rx{n) is a product of r over all nodes of the Young diagram of the partition A where 
argument of r is defined by entries i,j of a node [25]. The value oi j — i is zero on the main 
diagonal; the value j — i is called the content of the node. For zero partition one puts ro = 1. 
It was shown [27] that 

T,(n,t,t*)=^r;,(n)sA(t)sA(t*) (8) 
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(where the sum is going over all partitions including zero one) is a TL and a KP tau function, 
which we call the tau function of hypergeometric type. The Schur functions SA(t),SA(t*) are 
defined with the help of 

oo oo 

sx{t) = det/iAi-i+j(t)i<i,j</(A) , exp ^ z'^tm = J2 ^^^kiS) (9) 

m=l fc=l 

Here /^^(t) is called the elementary Schur function, or, the same, the complete symmetric 
function, see [25]. We should say more about the case t = t(x"). Mark that SA(t(x")) = if 
the length of the partition / (A) exceeds n, therefore, in this case, the sum over A (8) is restricted 
by < n. Due to (7), it means that r(/c), /c < are not involved to the series (8) if t = t(x"). 

If x"^ = (xi, . . . ,a;„) are eigenvalues of a n by n matrix X, and the variables t and x" 
are related by (5), we write Sx{X) := Sx(t{x.'^)). It is suitable to introduce the notion of tau 
function of matrix argument by analogy with the hypergeometric function of matrix argument 
(see (18), (20) below): 

r (n, X, t*) := r (n, t(x"), t*) , r (n, X, Y) := r (n, t(x"), t*(y'^)) , (10) 

where Hirota-Miwa variables xi, . . . ,Xn are just eigenvalues of a matrix X. We shall use large 
letters for this matrix argument. As in the case of the hypergeometric functions the tau function 
by definition depends only on eigenvalues of the matrix. 



2 Useful formulae 

We shall exploit the following formulae for integrations of Schur functions over the unitary 
group [25]. Let dM be the normahzed Haar measure on the U{n) and let S^^x be the Kronecker 
symbol. Then 



/ s,(AUBU-^)d,U='-^^^^, (11) 

•Ju{n) Sx{In) 



sx{A)sx{B) 

Sx{In) 

[ s,{AU)sx{U-'B)d.U = '-^^5,,x (12) 

JUin) S\[ln] 



s_x{AB) 

lUin) — ■ - ■ sx 

Let too = (1) 0, 0, 0, . . .). For integration over complex matrices Z there are the formulae [25] 



and 



/ ^ sMZ)s>.{Z^B)e-T^^^^ n = '-^S„x (14) 



Throughout the text d'^Z — tt "'^ Ilij=i dUZ^jd'^Zij] at sequel we shall omit C"^. Then we need 
[25] ^ 

SA(too) = ;^, Hx = l[i>'i + X'j-i-j + l) , (15) 

where Hx is the hook product. The Pochhammer's symbol related to a partition A = (Ai, . . . , A^) 
is the following product of the Pochhammer's symbols {a)x — {a)xi{a — 1)^2 • • • {a — k + 1)^^., 
{a)x, = r(a + A,)/r(a). Let t(a) = (f , f , f , . . .). Then 

{a)x = Hxsx{t{a)) = '^^^ , {n)x = HxSx{In) = (16) 
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In addition we have a simple relation (bosonic (l.h.s) and fermionic (r.h.s) representation of 
the vacuum TL tau function, since we did not find it in literature we proved it in [28]): 



(17) 



m=l 



which is a generalized version of Cauchy-Littlewood identity [25]. 

Let r{k) — q^. One gets tau function used in [29] as a generating function for double Hurwitz 
numbers. 

Let r{k) = Y[i=i{k + a^) n|=i(^ + and let /(A) be the number of non-vanishing parts 

of A. Let us write down the formulae for the hypergeometric function of a matrix argument X 
[11],[27] (throughout the paper all sums include zero partition) 



/ai + M,. . . ,ap + M 
bi + M,---,bs + M 



s 



X] =T,(M,t(x"),t„) = 



liX)<r 



Examples: 



oFoiX) = eTr^ , ,Fo{a\X) = e^^=^ ^^^"^ = det(l - X)' 



(18) 



(19) 



Now let r{k) — ni'=i(^ + ai){k + n — M) ^ni=i(^ + h) ^- The hypergeometric functions of 
two matrix arguments X,Y is [11],[27] 



fai + M,...,ap + M 



P"' s 



X,Y\ =Tr{M,t{^^),t*{y-))^ 



j2 nLi sx{t{ak + M)) sx{X)sx{Y) 



. ni=iSx(t(h + M)) 

i(A)<n 



sa(4) 



(20) 



Here = (a;i, . . . , x^), = (t/i, . . . , y„) are eigenvalues of the matrices X, Y. Matrices X, Y 
are supposed to be diagonalizable (i.e. normal) matrices. 

Taking bi = —c + e where c > —n is an integer, we obtain [18] 



lim rjFr 



e-»0 



p-t q 



Oi, . . . , Op 

-c + e, &2, ■ ■ ■ , bq 



(ai)a ■ ■ ■ (ap)a 



^detX^'+^Fq 



X 



ai + n + c, . . . , ttp + n + c 
2n + c,b2 + n + c, . . . ,bg + n + c 

where a is the partition (n + c, . . . , n + c) with the length l{a) = n. Meanwhile 



lim 



(i\ , . . . , dp 
-c + e,b2r ■ ■ :bq 



'-c + e)^{b2)a---{bq)^ 



{ai)a ■ ■ ■ (ap)o- 



detX"+'=detr'^+'^„JP; 



p^ q 



ai + n + c, . . . , Qp + n + c 
2n + c,b2 + n + c, . . . ,bq + n + c 



X,Y 



In particular 



lim iFi 

e->0 



n 

-c + e 



X)i-c + e)^ = oFq (X) det = e^^ det 
4 



(21) 
(22) 

(23) 
(24) 

(25) 



lim oFi 



a, n 
-c + e 



X 



-c + e)a 



detX^+^Fo {a + n + c\X) = det X"+'=det(l - X)""" 



(26) 

These relations are particular cases of the limit of r.,.. where we put r to be singular: r is 
evaluated at the vicinity of zero as r(e) = ^ + 0(e). Below c > —n is an integer and cr is a 
partition {n + c, . . . , n + c) of the length /(cr) = n. Then we have [18] 



r,(-c + e,t(x"),to,) detX"+'^ w nN , n -/^ ,m ^ 

lim ^ — ^ = —r-. Mn, t , too , nk) = r{k) , 

e-*o j-a{-c+e) [nja- k + u + c 



lim 

e->0 



T,(-c + e,t(x"),t*(y")) 



det X"+Met y"+V^(n, t(x''), t*(y'^)) 



(27) 



(28) 



r^(-c + e) 

Let us mark also the determinant representations of (18) and of (20) which are of 
importance in applications of the method of orthogonal polynomials to the new matrix models, 
which we shall consider below. 

The determinant representations have forms (the reader may find details in [16], [27]): 

det xrV(M-/c + l,t(x,),t*) I lk=i 



r,(M,t(x'^),t*) = 



A(x") 



(29) 



(which in particular yields the determinant representation of of (18)) and 

r,(M,t(x ),t (y )) = c,(M) A(x^)A(y^) ' 

(which yields the determinant representation oi pTg, see (35) below), where 

A(x")= n (^^ - ^,) = det(xr'=) , 

l<i<j<n 

n-1 

ci = 1 , Cn{M) = n {r{M -n + fc))'"" , n>l, 

k=l 

Tr{M -n+1, Xi, Vj) = 1 + r(M - n + l)xiyj + r(M - n + l)r(M - n + 2)x,2y| + 

where c~^(M)A(x")A(y") is supposed to be non-vanishing; otherwise the formula (30) should 
be modified, see [18]. Most often we need the cases M — Q,n. Examples of (29), (30): 

1 , , „ f ai — k + 1, . . . , a-p — k + 1 



(30) 

(31) 

(32) 
(33) 



/ ai, . . . , Op 



X 



^ f ai + n,...,ap + n 
^ \bi + n, - ■ ■ ,bs + n 



A(x") 

x,y] 



' " \h-k + l,---,hs-k + l 



Xi 



n 

i,k=l 1 



Cn{n) 



A(x")A(y^ 



■ det „F< 



ai + 1, . . . , Op + 1 



s 



fel + 1, 



(34) 



=1 (35) 



If r(0) = there is the following formula for the open Toda lattice (for instance see [42], [16]): 



/ Qi+k-2 \ 

Tr{n,t,t*) =Cndet [ a,^_ic,,,fc_i rr(l,t,t*) I 1^^^ , 



dt\-^dtl 
cn = n (r(fc)) 

k=l 



k—n 



(36) 
(37) 



At last let us notice that Tr(n, t,t*) enjoys the following symmetry relations: 

rr(n, t, t*) = Tj.(n, t*, t) , rr(n + m, t, t*) = Tf (n, t*, t) , (38) 

Tr{n,t^,t*) = Tri{n,t{a),t*) = Tr^{n,In,t*) , (39) 

where is n by n unit matrix, and where f{k) = r{k + m), ri{k) = {k + a)~^r{k), r2{k) — 
{k + ri)~^r{k), the relations (39) are derived from (16), for details see [27]. 



3 Tau functions of a matrix argument and angle integra- 
tion. The integration of tau functions over complex 
matrices 

Theorem 1 LetXi andyi, where i = 1, . . . ,n, are eigenvalues of matrices X andY respectively. 
If for n > 1 

n-l 

cn{n) = n ir{k))'-- + , (40) 

k=\ 

then there is the following generalization of HCIZ integral 

L) ""^^^^ ^ ^ ^"("^ A(x-) A(y») ' 

where A(x"), A(y") are the Vandermond determinants. 

The proof follows directly from (8), (11), (30). In particular we have (see (18), (20)) 



^ / fli, . . . , flp 



XUYU+ \dM^ 



Ol, . . . , Clp 



X,r (42) 



Let us consider the simplest examples of (42), which, at the same time, are the solvable 
cases of 



Similar series for matrix integrals were considered in the papers [5], [6] without connections with 
classical integrable systems. 

(1) HCIZ integral, [20],[21],[1],[30],[3], which we obtain choosing p = s = 

I Tr(xuYU+).rr V- Sx{X)sx{Y) v\ det(e^ 

/ e >^ M*C/= 2^ — —— ^ oJ^o{X ,Y) ^ Cn-r-—ri—r , (44) 

n—1 n— 1 / 1 \ k—n -i 

ci = 1 , c„ = c„(n) = n (r(^))'-" = n (^) ' ^ > 1 ' ^(^) = ^ (45) 

(2) Choosing p = s + l = lwe get 



/ det(l-XC/yC/+)~"d.C/= Y: {a)x'4P^-rr{n,X,Y) (46) 

Ju(n) V y ^^^^^^^^ i^ASA(4) 

- («, ^ , 1^) - A(x")A(y") ' (^^^ 

cn = c„(n) = n (K^))'"" , n>l, r{k) = (48) 
k=i 

Now let us consider integrals over complex matrices. Different integrals over complex ma- 
trices were considered in the papers of I.Kostov, see [31], [32] as an example. 
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Theorem 2 Under the conditions of the Theorem 1 we have 



j r. (n, t^, XZYZ^) e'Tr^^^d^Z ^ r. (n, X, Y) = ^^(^^ ^01(^1^.^))^,,...^^ ^^^^ 



where d'^Z = tt 11^7=1 d^Zijd'isZij. In particular we have 



jpFs 



ai, . . . , Op 



XZYZA e-Trzz^d-'Z = ,^,Ts r!"'' ' ' ' 



X,y] (50) 



The proof foUows from (8), (13), (30), (18), (20). Examples of solvable models provided by (50): 

(1) p = s = 0, then 

J,EZ.Jr(xzyz^)^-Trzz^^2z= ^ ..(X).,(F) = [1(1 - x.y,)"^ = i^o F) (51) 

\,l{\)<n i,j 

(2) p = s + l = 1: 

|eS:-^Tr(xzrz+)-^-'IVzz+^2^^ I _ ]^XZY Z^Y e-^^''^^ d'Z (52) 

= ^ (a)A5A(X)5A(r) = 2^o(«,^l^,n , (53) 

A,i(A)<n 

which in general is a divergent series (until —a is a non- negative integer). 

Next theorem is a generalization of Gross- Witten integral: 
Theorem 3 Let Zi,i — 1, . . . ,n be eigenvalues of the matrix XY 

j^^ ^ Tr {n, t, XU) Tf (n, In) dM = Trf {u, t, XY) (54) 

det (zr\r{n -k + l,t, t*(^i))) . 



A(z") 

where A(z") is the Vandermond determinant. Also 



(55) 



j^^ ^ Tr (n, t, XU) Tf (n, U-^X-\ t*) dM = {n, t, t*) , (56) 

where ri is the following step function: 

ri{k) = r{k)f{k) , k>0 , r{k) = , A; < (57) 

For the proof we use (12), (29). The simplest example of (56), X = l,r = f = lis the 
well-known model of unitary matrices, see [33] about this model. We obtain 

/ eS.*^Trt/^+S,t|Trf/-«^^^^ ^ 5,(t)5,(t*) = r.(n, t, t*) (58) 

X,l{X)<n 

In (58) r is the following step function: r{k) = 1, k > ; r{k) = 0, /c < 0. 

Due to the restriction l[\) < n it is not the l.h.s of Cauchy-Littlewood formula (17). Let 
us mark that the right hand side of (58) is the subject of the so-called Gessel theorem if one 



take t = t(x'"), t* = t*(y™), n < m. This case was considered in [4]. There exist different 
determinant representations of (58), which are due to (29), a modification of (30) and due to 
(36), see [18]. 

An example of (54) is 



r e^m^-^'^^"^^"^ det {l - U-'y)'' dM ^ E (b) 

X,l{X)<n 



sxit)s^{XY) 



(59) 



where Zi,i — 1, . . . ,n are eigenvalues of the matrix XY . Here hk{t) is the elementary Schur 
function, see (9). 

Consequence 1 Let A, X, Y are n by n normal matrices. For integer m > consider the 
partition a — (m, . . . , m) o/ the length l{a) — n. We have 

[ Tr {n, A, XU) Tf (n,U-% In) det U-'^dM = det A"" det X"'ra{n)Tr,{n, A, XY) , (61) 

Ju{n) ^ ' 

where ri{k) = r{k + m)r{k) . Also 

[ Tr {n, A, XU)rr(n,U-% In) det U"" dM ^detY"'r^{n)rr,{n, A, XY) , (62) 

Ju{n) ^ ' 

where ri{k) = r{k)f{k + m). In particular 



I 

JU{n) 



det IJ-'^nT. 



P"' s 



{n),X\U{h)c 



f 

JU(n 



fll, . . . , dp 

det/l™detXVp-^.+.-+i 

Cll , . . . , (Xp 



A,XU]pF, 



tti, . . . , Ojp 



U-^Y dJJ 



bi,---,bs 

ai + m, . . . , Qp + m, di, . . . , dp + m 
+ m,bi + m, . . . ,bs + m,bi, ■ ■ ■ ,bs 



A,XY 



IU{n) 
Ifi=l{ai)a 

{n).Uliik)a 



det V^pTs 



bi,---,b. 
det Y^pj^pJ- s+s+i 



A,XU\pFs 



(i\ , . . . , dp 



U-'Y dM 



ai, . . . , dp, ai + m, . . . , dp + m 
n + m,bi, . . . ,bs,bi + m, ■ ■ ■ ,bs + m 



A,XY 



(63) 



(64) 
(65) 

(66) 



Let us notice that putting A = in pJ^s{' " ■ l^j X) one gets pFs{- ■ ■ \X) with the same sets of 
indices {aj}, {bi}. Choosing o-Po of (19) we obtain the simplest example of (63)-(66): 



/ dett/^-e^(^^+^-^^)4C/ = (^"'^^2|oF,(n + m,XF) 
Ju{n) {n)a [ det Y"^ J 



(67) 



For the determinant representation see (29). At first this integral was evaluated in [34] via 
methods of [35]. The matrix integral in case m = was considered in [36] by a different 
method (by the method of orthogonal polynomials), and links with the KP equation and with 
the so-called generalized Kontsevich model were established. For m = and X = Y = —^In 
the integral is used for the study of two-dimensional QCD [22] and called Gross- Witten one 
plaquette model (about this integral see also [37]): 



Ju{n) 



(68) 



In n — >■ oo limit the Gross- Wittcn model enjoys a third order phase transition 
The other example of (63)- (66) is (cr is the partition (m, . . . , m), /(a) = n) 



(69) 

The determinant representation of the Gauss hypergeometric function see (29). Let us mark 

that 2-P'i(a., i>\ c\X) with integer a, 6, c solves Painleve V equation [38]. 

Let U (and V") be by (and respectively n by n) unitary matrices, N > n, and A, X 
(and B, Y) are n by iV (and respectively N by n) rectangle matrices . One combines results 
of Theorems 3 and 1 to get 

f f xuYV++Tr vau+b ^^^^^y ^ ^j,^^ri\BX, YA) (70) 

Juin) JU{N) 

At first the similar result was obtained in [34] (sec also [39]), where the answer was given as 
a determinant of Bcssel functions, by (35) it is the same as qJ-'i of (70). Other examples of 
combining of results of Theorems 1,3: 

/ / det(l - XVYV+r^'e^ ^^^^""dMd.V ^ iJ^i BX,Y a] (71) 

Ju(n) JU(N) \n J 



XY 



[ [ det(l - XUYV^)-" det(l - VXU^Y)-^dMd,V = aJ^i ( 

Ju{n) JU{N) \ n 



a, b 
n 



BX, YA (72) 



Theorem 4 Let Zi,i — 1, . . . ,n be eigenvalues of the matrix XY 

j Tr {n, t, XZ) Tf (n, Z+Y, too) e" '^^^^^d^Z (73) 



detf^;," ''Trf{n-k + l,t,t*{zi^ 
Tr, in, t, XY) = ^ '-^^^^^ (74) 



Also 



J Tr [n, t, XZ) Tf (n, Z+X-\ t*) e" ^^^^^d^Z = Tr, (n, t, t*) , (75) 
where ri is the following step function: 

r^{k) = kr{k)r{k) , k>0 , r{k) = , k<0 (76) 

The proof follows from (29). Examples: 

I eE™=i t,^Tnxzr^Tr{z+Y)^-Trzz+^2 ^ ^ ^Y.Z=i t-TV (77) 

/eS:^i*-Trz"^eS:^i*^'I^(^")'"e-^^^"d2Z = Y.{n)xSx{.t)s,{t*) (78) 

A 

Let us notice that series in the r.h.s. of (78) coincide with the series for hermitian-antihermitian 
matrix model and with the series for the model of normal matrices [17], which in n — > 00 limit 
describes the interface dynamics of a spot of water inside oil [40]. 
Applying (27) and (28) we get 
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Consequence 2 Let A, X, Y are n by n normal matrices. For integers m > we have 

j Tr (n, A, XZ) Tf (n, Z+Y, tj) det(Z+)"*e- '^'^^^^cfZ (79) 

= {n)a r^{n) det A"" det X"" r,, (n, A, XY) , (80) 

where the partition a = (m, . . . ,m) and are of the length l{a) = n. In (80) ri{k) = -^^^ r{k + 
m)f{k). Also 

j Tr {n, A, XZ) Tf (n, Z+Y, tj) det Z^'e' '^^^^^d'^Z = r^{n) det Y"" Tr, {n, A, XY) , (81) 

where the partition a — {m,...,m) and are of the length I {a) = n. In (81) ri{k) = r{k)f{k+m). 
In particular 



I*/ c 



Qj\^ . . . , Qj'p 



bi, ■■■,bs 



p-^ s 



Z+Yj det(Z+)-e-^^^^^ci2^ 
ai + m, Op + m, 5,1 ..... ttp 



ntl(^ det A- det X- , ^ ^ 

U.i=i{bi)a \bi + m,...,bs + m,bi,---,b 



A.XY 



I 



/ ai , . . . , Op 



Z+y\ detZ^e-^^^U^Z 



nLi(a»)^ det y"" -JF - ( «i' • • • ' «i + • • • ' «p 
Ut=i(bi)a \bi,...,bs,bi + m,---,bs + m 



A,XY 



For instance 



= detX^e^^^^ 



/det(z+)"eTr^^^ 
j det (z+)™det(l - Z+y)-«e'^^^-'^^^^d2Z = detX-det(l - XF)"" , 



(82) 

(83) 
(84) 
(85) 

(86) 
(87) 



j det det(l - XZ)-"" det(l - Z+Y)''e-^'^^^'' d^ Z = {a)^ det X"^2i^o(-A:, a|Xr) (88) 

Let Zi (and Z2) be N hy N (and respectively n by n) unitary matrices, N > n, and A,X 
(and _B, F) are n by (and respectively N by n) rectangle matrices. Theorems 4,2 yields 
the following examples: 

f f xziyz++IV z.Az+Bg-IVziZ+-IVz.z+^2^^^2^^ ^ o^o(SX, YA) (89) 

/ [ det(l - XZiyZ2+)-"eTi- ^2^^i+^e-Tr^i^i+-Tr22Z+^2^^^2^^ ^ ^qq^ 

4 Kontsevich-type integrals 

Lemma If /i, ...,/„ is a set of functions of one variable, and g{x."') is an antisymmetric function 
of xi, . . . , Xn, then 

J---J ^(x") det ifjixi)) \l^^,dx, ■■■dxr,^n\ j ■■■ j (^(x") f{ f^{xi)dx^ ■ ■ ■ dx„ (91) 



(below we put g(x^) = A(x")) 

Let X be n by n Hermitian matrix with eigenvalues Xi, . . . , x„, let and F be n by n matrix 
with eigenvalues yi, . . . , yn- Let dX = Hj^k d^Xi^d^Xi^ n"=i dX^^ = (i^[/A(x")^ 0^=1 dxi. The 
determinant representation (29), the lemma and theorems 1,2 result in 
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Theorem 5 The following integrals over Hermitian matrix X and over complex matrix Z 

J Tr(n, t, X)Tf(n, XY, QdX = J Tr(n, t, X)Tf(n, ZXZ+Y, too)e~^^^^^ d^ZdX = (92) 
is equal to 

n\C 



A(y-) 



— J det (a;^-Vr(n - A; + 1, t(xi), t*)Tf (1, Xi, Vi)).^^^ n ' ' ' ^^'^ ^^^^ 
det ( [ x'^~^Tr{n - A; + 1, t(x), t*)Tf.(l, x, yi)dx] (94) 

\J / i,k=l,...,n 



n\C 
~ A(y-) 

For the proof we use (41) and the Lemma (91) and (29). One obtains the Kontsevich integral 
|gt3Trx3g-Trxy^^ [41] if he choose r = 1, = and f{k) = 1/k. 
Also we have 

Theorem 6 Let ai, . . . ,an (yi, . . . , yn) are eigenvalues of matrix A (of matrix Y ) 

j Tr{n, A, X)Tf{n, XY, QdX = -^^-^ det (^j r^(l, a^-, a;)r^=(l, x, yi)dx^ ^ ^ (95) 

5 New multi-matrix models which can be solved by the 
method of orthogonal polynomials and the Schur func- 
tion expansion 

Here we consider in short new solvable multi-matrix integrals which we obtain with the help of 
the previous consideration, details will appear in the separate paper. 
Let us remind that the multi-matrix model of Hermitian matrices 

7 ^ y gTryi(Mi)+-+Try;v(Miv)gTrMiM2gTrM2M3 . . . ^TrM^v-iM^v^j^^^ . . . ^ (-gg^ 
where Mi, . . . , are Hermitian n by n matrices and 

oo 

V,iM,) = J2€^Mjr, k = l,...,N (97) 

m=l 

was shown to be multi-component KP tau function and solved by the method of orthogonal 
polynomials. 

Now we consider more general integral over matrices Mi, . . . , Mjv 

7 = y e'^^i^^i)+-+'^^~(^^)i?i(MiM2)i?2(M2M3) • • • RN-i(MN-iMN)dMi ■ ■ ■ dM^ (98) 

If we choose the interaction term as hypergeometric function of matrix argument (see (10)) 

Rk{MkMk+i) = Tr^{n, t, MfeMfe+i) , t = (1, 0, 0, 0, . . .) (99) 

then due to (41) it is possible to perform the angle integration over each Uk{n) where Mk — 
UkX^^^U^, X^^) = diag{xf\ . . . , x^f)). Really, (41) is a tau function which has a determi- 
nant representation (30). Finely we obtain the following integral over eigenvalues xf^"^ (i — 
l,...,n;k^l,...,N): 

n N 

I -of np.(^r^^fV2(xrxf ) • • •PAr-i(xr^^4^)) n e^t.^::-.^'-'^4''rd^p , (loo) 

•' i=l k=l 
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where 

p,(xf)xf^^)) =r,,(l,t,xfM'^^)) , t= (1,0,0,0,...) (101) 

(In the case of the famihar multi-matrix model (96) each pk{xy), k = 1, . . . , N is e~^'^). 

The integral (100) may be evaluated by the method of orthogonal polynomials {Pmif^ -.x)}, 
{Km{t^, x)} which depend on the collection of times t^ = (t^^^, t^^^, . . . , t^^^) 

j pm{t^, x)Trn{t^, vMt^, X, y)dxdx = (5„^e^»(*'') , n, m = 0, 1, 2, . . . (102) 

Each of polynomials is of a formp^(t^, x) = T.n<mPnm{t^)x'\ -Kmif^ , x) = J2n<m 7lnm{t^)x'^, 
and the weight function is 

. N-l 

uit'', X, y) = e^i(^)+^-(^) / pfe(xx(2))p2(x(')x(^)) • • • pAr-i(x(^-^V) 11 e^^^^^'^^dx^ (103) 

k=2 

We obtain 

7 = c n e'^'^*''^ (104) 

fe=i 

Let us notice that the orthogonal polynomials are Bakcr-Akhiezer functions of TL hierarchy. 
Gauss-Zakharov-Shabat factorization problem for TL equation [24] directly results from the 
relation (102) 

i^+(i^) = i^_(i^)e'^G(t^) , * = diag{ct>i{t^), Mt""), ■ ■ ■) , (105) 

where 

{K+)nm = (7r'~^)nm , {K-)nm = Pnm , Gnm = j X^'y'^Ujit^ , X, y)dxdx (106) 

(as it was in the case of two- matrix model [42]). 

The Schur function expansion for these matrix matrix models is given by (6), where K^^ 
are expressed via products of skew Schur functions, will be presented in the forthcoming paper. 
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